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We analyze the universality of the bosonization rules in non-relativistic fermionic systems 
in (2 + We show that, in the case of linear fermionic dispersion relations, a general 
fermionic theory can be mapped into a gauge theory in such a way that the fermionic 
density maps into a magnetic flux and the fermionic current maps into a transverse 
electric field. These are universal rules in the sense that they remain valid whatever the 
interaction considered. We also show that these rules are universal in the case of non- 
linear dispersion relations provided we consider only density-density interactions. We 
apply the functional bosonization formalism to a non-relativistic and non-local massive 
Thirring-like model and evaluate the spectrum of collective excitations in several limits. 
In the large mass limit, we are able to exactly calculate this spectrum for arbitrary 
density— density and current-current interactions. We also analyze the massless case 
and show that it has no collective excitations for any density-density potential in the 
Gaussian approximation. Moreover, the presence of current interactions may induce a 
gapless mode with a linear dispersion relation. 

1. Introduction 
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2 Functional Bosonization of Non-Relativistic Fermions 



The method of bosonization has proven to be a very useful technique in under- 
standing the non-perturbative regime of strongly correlated fermionic systems in 
one spatial dimension. Its analysis has shown that the bosonization rules are uni- 
versal, in the sense that they are model independent. Another important property 
is that the technique displays a powerful duality mechanism, which is evident in 
the connection of the non-perturbative (resp. perturbative) regime of the (1 -I- 1)- 
dimensional massive Thirring model to the perturbative (resp. non-perturbative) 
one of the corresponding bosonized sine-Gordon model EJ. 

In recent years, many efforts have been done in order to generalize bosonization 
to higher dimensions, with special emphasis on (2 -t- l)d fermionic systems, due to 
their relevance for low-dimensional condensed-matter physics as, for instance, the 
quantum Hall effect and high- Tc superconductivity. 

In the context of condensed-matter systems the first attemps to understand 
the bosonization of fiiictuations around the Fermi surfgice in higher dimensions was 
made by A. Luther □ followed by F. D. M. Haldane 0, who considered the two- 
dimensional system as a collection of infinite one-dimensional systems parametrized 
by an angle. Also, a functional bosonization method has been applied, essentially, 
to tha evaluation of single-particle Green's functions of interacting fermionic sys- 
tems □. The fermionic interaction is decoupled by making a Hubbard-Stratonovich 
transformation using a bosonic auxiliary field 4>{x). The Green's function is cal- 
culated in a fixed (j){x) background and then is averaged using an effective action 
that depends only on the Hubbard-Stratonovich field. This approach is suitable 
to implement several approximation schemesa and also to systematically analyze 
the case of fermionic systems with non-linear dispersion relationaj. Bosonization of 
Fermi liquids inhjgher dimensions has been also considered by A. H. Castro Neto 
and E. Fradkin Si. 

As for quantum field theory, two different approaches have been proposed to 
bosonize systems in higher dimensions, namely, the canonical and the functional 
methods. The canonical method was the one employed to obtain J;he first results 
for the bosonization of free massless fermions in two dimensions u. In this case, 
bosonization is achieved by introducing a vector gauge field, for which the bosonized 
action takes the form of a non-local Maxwell-Chern-Simons model. 

Functional bosonization, as a calculational tool, is undoubtedly a clear contri- 
bution to the understanding of strongly correlated electrons in higher dimensions. 
However, it is known in (1 -|- l)d systems that there are general properties that do 
not explicitly depend on the details of the one-particle Green's function, but do 
depend on global features such as symmetries and the universality of the fermion- 
boson mapping. For instance, it has been shown recently E3 that the conductance 
of finite incommensurate Peierls-Frohlich systems is universal due to a global chiral 
anomaly and to the universality of the (1 + l)d bosonization rules. For this rea- 
son, it is relevant to ask oneself whether there exist, in higher dimensions, similar 
universal rules. 
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A first step in this direction has been taken by E. Fradkin and F. A. Schapos- 
nik lii, who have analyzed the bosonization of a massive Thirring model in 2 + 1 
dimensions. They have shown that in the infinite-mass limit the fermionic current 
maps into a topological current, and the bosonized action is a local Maxwell-Chern- 
Simons one. Since then, a great deal of work has been done towards the understand- 
ing of the bosonization structure in 2 -I- 1 dimensions. While the canonical method 
is useful to investigate the mapping among the fundamental fields of the fermionic 
and the bosonic models, the functional method M0E3O't£l is appropriate to es- 
tablish the general framework for the higher-dimensional bosonization of current 
correlation functions. 

In a recent paper 0, it has been shown for relativistic systems that the mapping 
between a fermionic current and a topological bosonic one is, in fact, universal. This 
means that this mapping may be performed for any kind of current interaction. In 
the present article we extend these results to the case of non-relativistic interacting 
fermions. We show that this universality remains valid in the case of non-relativistic 
fermions with a linearized dispersion relation. In the case of a non-linear disper- 
sion relation, the boson-fermion mapping is universal, provided we only consider 
density-density interactions. As an example, we apply this technique to the case of a 
non-relativistic Thirring- like model with arbitrary density and current interactions. 
We compute the bosonized action in the Gaussian approximation and evaluate the 
spectrum of collective excitations of the model for arbitrary potentials. We find 
that in the large mass limit the system has collective excitations whose dispersion 
relation has a gap depending on the particular relation between the potentials. In 
the masless case, we show that the model has no collective excitations when consid- 
ering only density interactions. However, a gapless collective mode appears when 
current-current interactions are taken into account. 

The present article in organized as follows. In Section 2 we analyze the case 
of a linearized fermionic dispersion relation. In §2.1 we define a non-relativistic 
Thirring-like model and in §2.2 we review the functional bosonization technique. 
The Gaussian approximation for the Thirring-like model is developed in §2.3 and in 
§2.4 the two-particle correlation functions are evaluated. The spectrum of collective 
excitations for particular limits of the model are described in Section 3. Finally, 
in Section 4 we discuss the universality of the bosonization rules for interacting 
fermions with a non-linear dispersion relation and in Section 5 we present our 
conclusions. 

2. A non-relativistic Thirring-like model in (2 + l)d 



2.1. The model 

Let us present a non-relativistic (and non-local) Thirring-like model in two spa- 
tial dimensions. The non-relativistic character of the model is displayed in the 
kinetic part of the action as well as in the interaction terms. In order to develop 
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the functional bosonization technique presented in the following section, we find the 
Lagrangian formalism more suitable. For this reason, we present the action of the 
non-relativistic Thirring-like model using the imaginary-time (r) formalism at zero 
temperature (Euclidean metric) in the following form? 

S = So + S'int, (1) 

with 



So = j dr j d^x V (7/.9^ + m) V (2) 

and 



Sint = J dTidT2d^xd^y [Vi^o)ixi v) Joi^) Joiv) + V{x,y)Ji{x)Ji{y)] , (3) 
where the fermion field ip is written as 

The interaction piece of the action has been written in terms of currents defined 
as 

where we have adopted the representation 

70 = -10-3, 71 = viai, 72 = W20'2, (5) 

in which tri, (72 and are the Pauli matrices, vi and V2 could be identified with the 
components of some "Fermi velocity" , while Vo{x, y) and V{x, y) are symmetric bilo- 
cal arbitrary potentials describing density-density interactions and current-current 
interactions, respectively. The density-density potential represents a two-body scat- 
tering interaction and may take the form, for instance, of the usual non-relativistic 
Coulomb interaction. Although in a non-relativistic system the current -current po- 
tential is much weeker than the density-density one, we have introduced it here for 
two reasons. First, we have considered the potential V{x^ y) in order to turn the 
model more general, so that we could cover a greater set of situations by choosing, af- 
ter the calculations, different potentials (note that in the case oIVq = V = g^5{x~y) 
we recover the usual relativistic local Thirring model). Second, in the last few years 



several models of fermions coupled w 
tion of strongly correlated electrons 



s 



1 gauge fields were considered in the descrip- 
In particular, the Chern-Simons dynamics 



"In what follows we will use the following notation: Greek indices run from zero to two while Latin 
indices may take the values one and two. 
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for the gauge field, incorporated in the Chern-Simons-Landau-Ginzburg theory for 
the fractional quantum Hall effect Ea, induces long-distance current-current inter- 
actions iipon integration on the transverse gauge fields, which could be modeled by 
V{x,y) B 

It is not difficult to show that the kinetic term of the action can be formally 
manipulated almost in the same way that a relativistc one. In fact, the algebra of 
the above defined 7 matrices, 

{it^nA ^'i.Qt.i,., (6) 

induces a metric tensor g^u — diag(l, wi, 112). With this metric, the square of a 
three- vector k = (w, fci, fc2) reads — k^k^ = lj^ -l-w^fej Thus, we can work 

with a general Euclidean metric g^i, and perform, at the end of the calculations, an 
analytic continuationc3 of the metric coefficients vi and V2 ■ It is interesting to note 
that in the relativistic case of = = 1 this procedure is nothing but the rigorous 
definition of a Wick rotation. Moreover, we shall find interesting topological results 
that do not depend on metric, so that without lost of generality we can work with 
vi = V2 = \ and make, at the end of the calculation, the correponding replacements. 

The model presented in this section was extensively analyz ed in (1 4- l)d in the 
context of the study of many-body physics of quantum wires EilHSa. It is the 
purpose of the present paper to extend its analysis to 2 + 1 dimensions . 

2.2. The bosonization technique 

In this section we present a functional bosonization technique in 2-1-1 dimensions 
which shall be applied to the particular case of the model described in the previous 
subsection. This technique is suitable for the calculation of current correlation 
functions, as it essentially maps a conserved fermionic current into a topological 
bosonic one. In 1 -I- 1 dimensions the bosonic field is a scalar one, provided the 
currents map as 

= V'TmV' — ' = ^t^'^d^'P- (J) 

Note that is a conserved current [d^J^ = 0) due to the global U{1) invariance of 
the theory; in other words, it is conserved due to the field equations of motion. On 
the other hand, the conservation of the bosonic version of the current {d^J^ — 0) 
is automatic and does not depend on the equations of motion. In this sense we say 
that the current is bosonized into a topological one. 

If one wishes to apply the bosonization technique in 2 + 1 dimensions in an 
analogous fashion, one has to remember the fact that the only way to build up a 
topological current in 2 -I- 1 dimensions is by employing a vector field, instead of 
a scalar one. Thus, we are looking for a bosonization scheme that maps fermionic 
currents into bosonic ones in the following way: 

= V'TmV' — ' = f^t^upd^Ap. (8) 
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Note that, again in this case, 9^ = foUows automatically, due to the totally 
antisymmetric Levi-Civita tensor e^i^p in (2 + l)d. It is interesting to remark that 
the fermionic current is invariant under a local U{1) gauge transformation of the 
fields given by 



— ' e*"(^V (9) 
^ _^ ^e-'"^^^; (10) 

conversely, the bosonic topological current also has this invariance, since, when the 
vector field transforms as 

— > Ap+idf,a{x), (11) 

then 

J'p = efj,^pdi,{Ap + idpa{x)) = Jp, (12) 

due to the antisymmetry of ep,vp- Therefore, we expect that, after the bosonization 
process, the bosonized theory to be a gauge theory. 

In the functional formalism, the current correlation functions are evaluated from 
the generating functional 

Z[s\ = j V^Vi;e-^^^'^^-'I'^''' '^^^ (13) 
by functionally differentiating equation (|l^ ) with respect to the source s^: 

Jiffl^. (14) 

The aim of this program is to find a bosonic representation to the generating func- 
tional of current correlation functions in the form 

Z[s\ = j 25^^e-^?bo3[A,]-»/rf^. s,j^^ ^^5^ 

so that 

By comparing eqs. ( [l^ ) and (|l6|) we see that the theories described by S[ip] and 
S'bosi^p] are equivalent in the sense that 

(Jj^(xi) . . . j;„(2:„))5 = {J^S^i) . . . J^^„(x„))5,„.. (17) 

In order to write a formal expression for the bosonized action and currents, we 
follow the path-integral approach of refs. E]E3ll3. We can write the generating 
functional for interacting massive fermions as 
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Z[s\ = j V^jV^/j e"/'^'" ^(^+n^+^M-S'[J]^ (^g) 

where we use the notation ^ — ^f^s^ and ^ = jf^d^; S^[J] is an arbitrary interaction 
action depending only on fermionic currents. The non-relativistic Thirring model 
presented above is a particular case corresponding to quadratic current interactions. 
In what follows we shall consider any type of current interactions (not necessarily 
quadratic) . The integration of the fermionic variables is a very difficult task due to 
the interaction part of the action which may contain non-quadratic fermionic terms. 
The usual trick in the case of local quartic interactions is to perform a Hubbard- 
Stratonovich transformation in order to obtain a quadratic fermionic action, at the 
cost of introducing a new auxiliary field. However, for an arbitrary interaction, 
such a transformation is not available. Here, we shall obtain a quadratic fermionic 
action for arbitrary interacting potentials, provided that they only depend on the 
fermionic current. The technique that we will use was presented for the first time 
in ref. 0. 

We begin by reexpressing the exponential of the fermionic interaction S^[J] in 
terms of its functional Fourier transform. 

Therefore, for the generating functional of eq.(|l^), we get 

In this way, we can interpret the fermionic interaction as an effective interacting 
fermionic system minimally coupled to a dynamical quantum vector field whose 
action is S[a]. Note that for a generic interaction S^[J] which is not quadratic 
in the current the computation of S[a] is, in general, not possible. However, 
we will see that our final result will be given only in terms of the known quantity 
^^[J], without relying on the explicit computation of its Fourier transform S[a]. In 
other words, all that will be needed is the reasonable assumption that the Fourier 
representation ( p^ ) does in fact exist. In particular, for a local density-density 
two-body interaction, S[a\ is quadratic and ( p^ ) coincides with the usual Hubbard- 
Stratonovich transformation. 

The next step towards bosonization is to decouple the external source from 
the fermion fields. This is accomplished by first making the following change of 
variables (with trivial Jacobian) 



^ — > V;e-*"(^\ (21) 
Vi(jVi} — > Vi}V^. (22) 
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The resulting generating functional Z[s] does not depend on a; therefore, we can 
integrate over a modifying only a global multiplicative factor that will not enter 
the correlation functions. We thus have 

Z[s] = J VijV^pVa^Va x e" / 'A(#+«(^+^+#")+'")V'-S[a] ^23) 

Let us now introduce the field bfj_ through the condition 

d^a = b^. (24) 

In order to integrate over &^ (instead of a) in the path-integral ( p3| ) we must impose 
that bf^ be a pure-gauge field. This can be done by inserting the delta functional 
H^fj.i^pfi'pM) the expression for Z[s], f^p — di^bp — Opb^, being the b field strength. 
Thus, 

Z[s] - J V^Vi)Va^Vbp 5{ep,pUp[b]) e~ ^{9+^{^+i+1>)+rny^-s[a\ ^ ^25) 
The source now decouples from the fermions by shifting the field through bp 

bfi — > b^- Sp- Gp, (26) 

so that we obtain 

Z[s] = [ V^Wi>VapVbp 5{ep,pUp[b - s - a]) e"/^'^ ^m^Hm)i,-s[a] _ (27) 



The fermionic variables now may be integrated out since their integral is quadratic 
and is decoupled from the source Sp. The resulting fermionic determinant may be 
written as an exponential factor, through the property IndetO = TrlnO, and we 
get 

Z[s] = j VapVbp S{ep,pUp[b - s - a]) eTrin(^+,R™)-S[a]^ (28) 

Exponentiating the delta functional by means of a Lagrange multiplier Ap, the 
generating functional becomes 

Z[s] ^ JvapVAp e-Ku[A]~^jd^x e,^,A,dAs,+a,)~S[a]^ ^29) 



where 



isolates the 6-dependent part of the generating functional. We now see that the 
exponential of the bosonizing free action Kb is obtained from the transverse Fourier 
transform of the exponential of the effective action Trln(^ + ip + m). 
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Fermionic 


Bosonic 


Fundamental Fields 




A, 


Conserved Current 




f) A 


Kinetic Term 


J cPx ip{J^ + m)il) 




Interaction Term 







Table 1. Bosonization rules for fermions with linearized dispersion relation in 2 + 1 dimensions 



Integrating now over and using again the representation (|l9|), we finally 
obtain 

Z[s] ^ J VAf, g-Swl^Ml-'/'i'^ s^e^^.d^A,^ ^3^^ 

where 

Sbos[Ap]=KB[A] + S'[edA]. (32) 

As anticipated, we see that our final result relies only on the known action , the 
unknown action S[a\ being only needed at the intermediate steps. We also see that 
the bosonized theory is a gauge theory. This statement is simple to verify. The 
interaction part of the action depends only on the current that is automatically 
gauge invariant as shown in (|l^. The gauge invariance of the kinetic term Kb 
can be understood from equation (^0|). If we make a gauge transformation A^^ — > 
Afi + df^a, the last term in the exponential of (|30| ) transforms as 



J d^x e^up {A^ + d^a) d^bp = J d^x SpiypAf^d^bp - J d^x aSpupdpdybp (33) 

where we have integrated by parts the last term. Of course, this term vanishes, due 
to the antisymmetry of Sp^p, from which results the gauge invariance of Kb- 

Comparing eqs. (^8|) and (|3l|), we can read off the bosonization rules for the 
fermionic action as well as for the fermionic current. We see that the bosonic field 
has the property of a gauge field, and that the bosonic current takes the form 

Jp = ^pupduAp \ J _ ^ ■ (34) 
I 'J I — '^ij-^j 

Therefore, the fermionic density bosonizes to a "magnetic" field B and the fermionic 
current to a transverse "electric" field Ei, both related to the ^^-field. These are 
the main results of this section and they are summarized in Table ^ 

Another interesting and important result is that the kinetic term and the inter- 
action term of the fermionic action are bosonized independently. That is to say, the 
bosonization process does not mix these terms. Last but not least, the equivalence 
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between the fermionic and the bosonic current is not only exact but it is also a 
solid universal result, in the sense that it remain the same in the interacting quan- 
tum theory. We have shown that we can bosonize any term involving the fermionic 
current by simply replacing by the bosonic topological current e^i,pdjjAp. 

The universality for such a large class of interaction terms shown in this sec- 
tion is strongly dependent on the fact that the kinetic theory has only first order 
derivatives. We shall come back to this point in §4.2. 

2.3. The Gaussian approximation 

As we have seen, the results depicted in table ^ are exact and universal. In other 
words, they do not depend on any particular detail of the fermionic interaction 
provided that it depends only on the fermionic vector bilinear. In this sense, we 
have exactly bosonized the current and a general class of interactions. However, 
the kinetic fermionic action is bosonized into /^^[A^] given by eq. (^0|). Thus, 
in order to obtain an explicit expression for this action we have to evaluate the 
transverse Fourier transform of the exponential of the fermionic determinant. Of 
course, this is a very difficult task and no exact result is known in two spatial 
dimensions. Nevertheless, expressions (|30|), ( |3l| ) and ( |3^ ) constitute the starting 
point for different approximation schemes not always available in the fermionic 
version. For example, the semiclassical approximation is very easily performed in a 
bosonic theory, while it is very hard to handle in the fermionic counterpart. 

A very interesting observation about the action Kb was made in ref. 0. It was 
shown that there is a non-linear and non-local field transformation that maps the 
bosonized action Kb into a pure local Chern-Simons form. The underlying topo- 
logical structure of the bosonized action is then made explicit and, in some sense, 
this transformation is a convenient form to distinguish the topological contents of 
a fermionic theory from the non-topological dynamical contributions. However, 
although this formal observation is a clear positive contribution towards the under- 
standing of the bosonized action, the practical applications of such a structure as a 
reliable computational tool remain undeveloped. 

Therefore, it is necessary to develop an approximation scheme for the evaluation 
of current correlation functions. In ref. 113, the limit in which the "mass" parameter 
m is taken to infinity was considered. This is equivalent to a low-energy expansion 
of the fermionic theory, since the energy fluctuations must be much smaller than 
the typical energy of the gap, 2m. In this limit, the bosonized kinetic term takes 
the form of a local Maxwell-Chern-Simons term: 

/i 1 
d^x-e^^pApd^Ap + —Fp^Ff,^ + ■■■. (35) 

This is a quadratic action, simplifying in this way the evaluation of current correla- 
tion functions. It is worth stressing that in this limit the bosonized action is exact, 
in the sense that the next quartic term is proportional to 1/m? . 

In ref. El, a more general approximation was discussed, the so-called Gaus- 
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sian approximation, in which one takes into account the fuU quadratic part of the 
fermionic determinant, which may be written as 



with 



Tr\n{^ + ip + m) = Tpcib) + Tpv{h) 



(36) 



Tpc{b) = I d'xF,,{b)F{~d^)F^,{b), 

Tpv{b) = [ £xb^G[^d^)e^,,xd,bx, 



(37) 
(38) 



where Tpc and Tpy come from the parity-conserving and parity-violating pieces of 
the vacuum-polarization tensor, respectively!!. 

A standard one-loop calculation for the Fourier transforms of the functions F 
and G in (H) yields 



F{e) 

G{k') 



Am? 



^ Am? ^ 



arcsin(l 



m Am, 
arcsm(l H — -— ) 



27r I k 



k^ 



Am^ 



(39) 



where here and in what follows we shall always denote momentum-space represen- 
tations by putting a tilde over the corresponding coordinate-space representation 
quantity. 

In order to evaluate Kb, it is necessary to calculate the functional transverse 
Fourier transform of (^). Since this expression is quadratic, it is a simple task to 
evaluate the functional integral provided we have properly fixed the gauge. We thus 
obtain 



j d^x I ^F^, Ci F^, - ^Af, C2 e^.A^, | , (40) 

where 

F 
G 

C2 = _d2F2 + G^- (42) 

In the limit to — * oo we recover expression (|35|). Although in this limit the 
quadratic approximation is exact, it is not obvious how to determine its range of 
validity for small to. However it was shown in refs. and that this quadratic 
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approximation is in fact the exact one for any "mass" , provided we use it to calculate 
two-point correlations functions. For higher order correlation functions it is not at 
all easy to justify this approximation for small to. 

Turning now to the interaction term of the model, its bosonization is performed 
by just replacing Jq ^-^d Ji — > \Ufj.vF^i, in (|^). We thus obtain 



S'[A^] = / cPxd'y {e,jF,,{x)V(o){x,y)ei^F,r^{y) + F,^{x)V(^){x,y)F,o{v)] . (43) 



Finally, from eqs. ( pO| ) and (|3|) we can write the bosonized action for the non- 
relativistic Thirring model in (2 + 1) dimensions in the Gaussian approximation 
as 



(fx(fy \ 2<^ijF,j(x)VQ{x,y)ei„,Fi„i{y) + FM{x)V{x,y)Fio{y) ^4) 



2.4. Current density correlation functions 

In order to calculate current-current correlation functions, we need to evaluate 
the generating functional Z[s] given by (|3l|) and (|32|). We shall use for S'bos the 



Gaussian approximation obtained in (44) 



Since we are dealing with a gauge theory, to functional integrate the generating 
functional we have to properly fix the gauge in (^4|). We gain more insight of the 
problem by choosing the Coulomb gauge V • ^ = and rcexpressing the gauge-fixed 
action in terms of Aq and the associated magnetic field B = V x A. 

We note that, since V • ^ = 0, we can always write the components of the field 
A as 

A, = ey9,0, (45) 

where (j){x) is an arbitrary scalar field. The relation between B and A implies, 
together with (^5|), that B is none other that the Laplacian of the field cj) (apart 
from a sign), so that the vector potential is expressed in terms of B using the 
Laplacian Green's function: 



A = -e,,^B. (46) 



In this way we can write the generating functional in the Coulomb gauge as a 
function of and B: 



Z[s] = I VAoVB f,'SUAo,B)~iJ d'^x soB+s,{e,,d,Ao-^B) ^ ^^^^ 
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with 



5, 



bos 



- / d^xd^y 



B{x) i^C,{l + ^^ + {Vo + V^^ I (a; - y)B{y) 

, (48) 



-Ao{x) (Ci + V)]{x- y)A^{y) - iAo{x)C2{x - y)B{y) 



where Ci and C2 are given by (|4^) and (^2|) and Vo{x — y) and — y) are the 
density-density and current-current interaction potentials, respectively. 

At this stage of the calculation, it is convenient to write the action in momentum 
space. Thus, by Fourier transforming (Esl), we obtain 



5b 



1 f dojdk ~, \ ~ / ~ \ ~ 

2i 27r ' 11 |fc|2 / I |fc|2/ I 



1 f diodk ~ 
'2 J (2 

diodk ~ 



(2") 



3 Ao{uo,k)C2B{-u,^k). 



Integrating out the Aq and B fields, we find 



where the vacuum polarization tensor 11^^ takes the form 
n^^^(a;, k) = Viu), k) ( 5^,^ - 



) +iQ{uj,k)ef,^pkp, 



(49) 



(50) 



(51) 



with 



and 



V{uj,k) = {u;^ + k^)(v + Ci'^A{uj,k) 



A{uj, k) 



V + C, 



Ci+V] +k^{Ci+Vc 



CI 



(52) 
(53) 

(54) 



The vacuum polarization tensor for the interacting theory II^i, has the correct 
symmetry properties. In fact, 

iipAu,k) = nt^{io,k) ^n,^{-co,-ky, (55) 

also, it must be transverse, as a consequence of gauge invariance: 
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kf^Uf^„ = n^^^fc^ = 0. (56) 

This condition is automatically satisfied due to the tensor structure of 11^^ (see 
eq. d^)). 

It is interesting to note that we can rewrite the parity-conserving part of the 
vacuum polarization tensor as 

n^^'^) = -^ — ^ ■ (57) 



This formula should be compared with the one-spatial-dimensional result M. We 
see that the effect of the kinematics in two space dimensions, given by Ci and C2, is 
twofold. The last term in the denominator (proportional to C2) is due to a parity- 
breaking effect and we shall show that it is responsible for opening gaps in the 
excitation spectrum (see Section 3). The parity-conserving effects renormalize the 
potentials as — and Vq^^ = . Note that this is a dynamical renormalization 

(since Ci depends on w), which changes the structure of the interactions. However, 
we shall see that, in the large in limit Ci is a constant and the net effect is just a 
coupling-constant renormalization. 

Taking into account this structure for H^^ , it is now simple to write an expression 
for the current expectation value as we have, by definition, 

W^^MW; (58, 

from eq. ( |50| ) , it is expressed in terms of II^i^ in the form 

(J^) ^ U^^s^ = n^oso + n^iSj. (59) 

Gauge invariance implies that the components of H^^ are not independent (see 
eq. ([5^)). In particular, we have 

n^o = --I^^.^k^. (60) 

Using this relation and defining an external electric field as E°^^ = Vsq — Oqs, we 
find 

{J,)=^-U,,E]-\ (61) 

UJ 

It is interesting to have an explicit equation for the density fluctuations (in the 
linear response approximation) induced by an arbitrary external electromagnetic 
field. We obtain from (|l|) and (|6l|) 
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' r{Lo, k) hE^' + - QiLO, k) e^jhE'^^' 



+ fc2 U! 

= -A{uj, k)(V + Ci)T (V • _ ^ ^^AO^, fc) ^ ^ ^cxt^ ^ (g2) 

where . .) denotes Fourier transform. 

Noting that the external classical field E°'^'^ satisfies Maxwell's equations, we 
finally obtain 

(p) = ~K{u:, k) [p-^\V + Ci) + C2B''^'] . (63) 

We see that a density-charge excitation in the system may arise as a consequence 
of external charges and/or external magnetic fields. The general relation between 
a magnetic field and density comes, in general, from the parity-breaking terms of 
the model. In the fermionic version this is seen as a gap in the spectrum and, in 
the bosonic picture, as a Chern-Sinions-like action. Equation (|6^ ) is a generaliza- 
tion of the Aharonov-Casher calculation El to interacting systems submitted to an 
arbitrary electromagnetic field. 

Similarly to (p), we can write a expression for (Ji): 

\ ^ J Uj{LJ^ + k^) 

In the case of a divergenceless external electric field, we find 

{Muj,k)) = a,j{LO,k)EY\uj,k), (65) 
where we have introduced a conductivity tensor 

for whose elements we have 



crj;a;(w, fc) = cr^^(^, fc) = '^(^' , (67) 

a^yiio.k) = QiLu,k). (68) 

In the same way we have calculated expressions for the current mean value, we 
can write expressions for the two-point current correlation function, as it takes the 
simple form 
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In particular, using (pT|), the density-density correlation function may be written 
as 

{pp)=Ii(,o = l?{V + Ci)K{uj,k). (70) 

The excitation spectrum is given by the analytic properties of A. In the next section 
we shall develop this point in detail. 

Similarly, the spatial components of the two-point current correlation function 
also take a simple form: 

(JiJi) = (F + (7i)(co2 + fc|)A(cc>,fc), (71) 
(J2J2) = {V + Ci){u;'' + kl)K{u;,k), (72) 
(J1J2) = - {(F + COfcifca + C'2t^} A(tj,fc). (73) 



3. Spectrum of collective excitations 



3.1. The large m limit 

We have seen from the density-density correlation function (eq. ([70|)) that the 
spectrum of collective excitations can be read off from the singularities of A(u;, k) 
in the complex u> plane. Let us analyze here the case of a large fermionic gap where 
the simplest case to be considered is the infinite-gap limit. Performing the m ^ 00 
limit in eq. (^9|), we obtain 

lim F{k) = 0, (74) 
lim G{k) = -L, (75) 

implying Ci = and C2 = (l/27r)'^. Replacing these values into eq. (jsj), we obtain 

A(Lu,k) = ^ . (76) 

I p.Vo I 2 ^ ^ 

In the infinite-gap limit in the absence of current-current interactions, A does not 
depend on uj, which implies the nonexistence of any kind of density excitation. 
Moreover, the conductivity of the system is (t^^ = (^yy = 0, while axy = 1/2 (in units 
of e^/h). This result is completely independent of the density-density potential Vq. 
The presence of a transverse current comes from the fact that the fermionic system 
is not invariant under a parity transformation due to the term proportional to 
m^ip. As a consequence, the system behaves as an insulator {axx — 0) whatever 
the potential Vq; although there is a transverse current, it does not dissipate energy 
since it is perpendicular to the applied electric field. It is important to stress that 
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this is the same behavior of free fermions with a large gap, as the effective action for 
the external electromagnetic field is a pure topological Chern-Simons one. Through 
this simple exercise, we thus verify that the same system but with an arbitrary two- 
body interaction (p-p) remains topological in the infinite-gap limit. 

Turning back to the case in which the current-current interaction is present, 
still in the infinite gap limit, we realize that it plays a crucial role. The potential V 
induces a pole in the density-density correlation function associated to a collective 
excitation with the dispersion relation 




(77) 



and residue 



Res(A) = — . (78) 

Consider, for example, the case of a local current-current interaction with coupling 
constant and a Coulomb density-density interaction Vq = /\k\. In this case, 
the collective excitations have a gap proportional to so that, the stronger the 

current interactions, the weaker the energy necessary to excite the collective modes. 
A direct consequence of the existence of this gap is that the conductivity is not 
affected by these excitations. It is simple to verify that also with a strong current 
interaction, a^x — and axy — 1/2. This topological property breaks down if the 
gap disappears as it happens with a potential Vq = l/|fc|" with a >2. 

A different situation emerges when the fermionic gap 2m is large but not infinity. 



We now keep a 1/ m-contribution from eq. (|39|) , so that A is given by eq. ( 54 ) with 
Ci — i2xi6TT^m ^^'^ ^2 — (l/^'"')'^- In these conditions, the espectrum of collective 
excitations contains a mode with dispersion relation 



±M±^P+ , (79) 



and residue 



Res(A) = — . (80) 

As expected, this case has a similar behavior to the previous one (m oo) with the 
only difference that the potentials Vq and Vi are renormalized by the constant Ci. 
In fact, the parity-conserving contribution of the dynamics of the free fermions acts 
as local density-density and current-current interactions with coupling constant Ci 
renormalizing in this way the "true" potentials. As a final remark, note that these 
results are exact for any interaction. In particular, eq. (^9h represents the exact 



18 Functional Bosonization of Non- Relativistic Fermions 



dispersion relation of the collective modes for any arbitrary static potentials Vo(fc) 
and V{k). 

3.2. The massless case 

In this section we analyse the m — ^ limit of the non-relativistic Thirring model. 
Making m — > in (BS) we find 



\miF{uj,k) = (81) 



implying that 



lim G{uj, k) = 0, (82) 

m— ^0 



limCi(tJ,fc) = (83) 



This result leads to 



limC2(w,fc) = 0. (84) 



A- . —. (85) 



Here, several comments are in order. First of all, expression ( pSj ) depends es- 
sentially orieq. ( p^ ) , which is not well established and it seems to be regularization 
dependent cJ. In ref. ts, we argue that if we use a regularization that does not 
introduce any additional symmetry breaking on the system, the correct result is 



that of eq. (82). Fortunately, although the value of A depends on the particular 



value of limm^o G, it was shown in ref. lj that the collective-excitation spectrum 
does not depend on it, being insensible to the parity-breaking term of the ferniionic 
determinant. 

To understand the excitation spectrum, let us simplify expression (85) by con- 
sidering a pure density-density interaction (V — Q). In this case, eq. (85) reduces 
to 

A = ^ . , ■ (86) 

This equation for A has no poles. Instead, it has a cut for k"^ < uj'^ whose disconti- 
nuity 

5+ A = -i lim I A(w + ie, k) - K{uj - ie, fc)| (87) 

is given by 
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where is the Heaviside function. This result shows that the model in the zero 
mass limit (in the Gaussian approximation) has no collective excitations. Instead, 
there is a continuum of modes that contribute to any physical observable. For 



instance, from eq. (63) it is possible to calculate the density fluctuations induced 
by an external charge. However, this density profile cannot be interpreted as a 
superposition of density waves with a definite dispersion relation 0. 

Let us consider now the effect of a current-current interaction on this system. 
For simplicity, let us consider Vq = and V — in (^5|). We immediately realize 
that we have two regimes of momenta. For fc^ < uj'^ we have a cut in the spectrum 
whose discontinuity is given by 



Moreover, for high momenta k'^ > uj'^, the current- current potential induces a 
gapless mode with linear dispersion relation uj — v\k\ with v = i_f_g2/iQ2 < 1- 

4. Non-relativistic spinless electrons in (2 + l)d 

4.1. The current generating functional 

In this section we consider a spinless non-relativistic fermionic system in two 
spatial dimensions with kinetic action given by 



(90) 



This action has a global U{1) symmetry and as a consequence it has conserved 
Noether currents. 



Jo = r^, (91) 

J, = V'*V,V, (92) 

in such a way that 9t Jq + V • J = 0. 

In principle, any functional of the currents Si{jQ,Ji), could be considered for 
the interaction part of the action, thus we can write quite generally 

S[r,i'] = So[r,i'] + SI[Jo,J^] (93) 

The idea is to develop the bosonization technique discussed is the last section to 
evaluate density-current correlation functions. The main point of this technique is 
to couple the system to an external electromagnetic field and use it as a source to 
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evaluate current expectation values. In the present case, this is not so simple to 
accomplish, due to the non- linear dispersion relation [uj = fc^/2m). More explicitly, 
let us couple the action ( p3| ) to a vector gauge field (sq, s). The result is 

j d^xdt ^* + so + -^^2^ - m| ^ + Si{Jo, J + sJo). (94) 

The main difference with respect to the linear case is that here the interaction 
action is modified due to the derivative operator that appears in the current. In 
other words, Ji is not a local gauge-invariant current and therefore it picks up a 
factor Si which couples to Jq. By functional differentiating ( |9^ ) with relation to 
(so, Si), we find 

^ ^^^^ 

SS ^ Jq SSr ^ 

^ = ^^+2^^' + ^-^- 

We see from ( p5| ) that in the limit ^ the potential sq may be interpreted as a 
source of Jq. However, from ( p6| ) it can be deduced that in the same limit Si only 
acts as a source of Ji provided 

SSijJo^ = 0, (97) 
oJi 

implying that, in order to use the functional bosonization technique developed in the 
present article for the evaluation of current correlation functions, we must consider 
only density-density interactions. 

4.2. The bosonization technique 

Consider the generating functional 

Z[s] ^ J VipV^* e-5['^*''^'«^'l, (98) 

where 

S[r,^, sp] = J d'xdt r + so + ''^^^^ - Aij ^ + Si[Jq\. (99) 
In order to integrate the fermions, we first rewrite the interaction action as 

SO that 

Z[s] ^ J V^ViP*V4> e-S['A*.'/'>V.0] (101) 
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Fermionic 


Bosonic 


Fields 




A, 


Density 




B 


Current 






Kinetic Term 


Jdxr{^^t + ^~^^}^p 


Kb{A^) 


Density Int. 




S'{B) 


Current Int. 




Not Universal 



Table 2. Bosonization rules for fermions with quadratic dispersion relation in 2 + 1 dimensions 



where now we have 

S = j d^xdti^* l^dt + So + 2r^^ - /^| V' + S[(j)\. (102) 

The decouphng of the source from the fermions and expressing the generating 
functional as a function of a gauge field follow from the same formal mathemat- 
ical steps that one takes from eq. ( p2| ) to eq. (^2|), obtaining 

Z[s\ = j VA^ g-Sboa^Ml-'/rf'^ soB+s,e.,E,^ ^^Q3^ 



S^os[A^] = Kb[A]+Si[B] (104) 



, Tr In ( idt +bo+ '"^otf* ti] +i f xe,,u pA,,d,jb„ 

Kb[A] = X)b^e \ 2m ) J ^^pg^ 



where 



and 



In this way, for a non-relativistic gas with density-density interactions, we find 
bosonization rules very similar to the ones encountered for the fermionic system with 
linear dispersion relation. In fact, the exponential of the bosonized free ferniion 
action is the transverse Fourier transform of the fermionic determinant and the 
bosonized interaction action has the same form as the fermionic interacting action 
when substituting the density p by the magnetic field B. Moreover, from the cou- 
pling to an external source s^, it is simple to read off the bosonization rules for the 
currents: 

= ^^.upd.Ap I t - p • (106) 

This means that we can map any fermionic theory in 2d into a gauge theory by 
associating charges to magnetic fluxes and currents to transverse electric fields. 
These results are summarized in Table 2. 
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Eq. (105) for Kb is an exact formal result and can be used to develop different 
kind of approximation schemes. As an example, let us consider here the Gaussian 
approximation. 

In the Gaussian approximation the fermionic determinant takes the form 



Trln ( i9t + 60 + ^ (V + - ) = 



n j cPx boix) - J <Px<Py bf,{x)ICf,^{x - y)b^{y), (107) 

where n is the mean density of electrons. The first (linear) term is a tadpole 
contribution and can be dropped out simply by measuring densities from n. The 
second term, the non-relativistic vacuum polarization tensor, gets the contribution 
from two types of diagrams, the bubble and the tadpole diagrams arising form 
the interaction b'^xjj*ip. It is simple to evaluate Kg, since the transverse Fourier 
transform is quadratic. Thus, after a proper gauge fixing is introduced, we can 
perform the functional integral, obtaining 



We have used the fact that the vacuum polarization tensor is symmetric in this 
approximation, which implies that the free system does not break parity. For this 
reason, a cross-term BEi is absent from Kb- However, it may happen that con- 
figurations with strong magnetic fields, that e xpl icitly break parity, could induce 
a Chern-Simons-type action. In refs. E3 and c3, this behavior of the functional 
determinant was deduced for long distances; the resulting action is non-quadratic, 
making the evaluation of the transverse Fourier transform a very difficult task. 
However, it is possible to perform it in the limit of a strong magnetic field, in which 
it is possible to project the system onto the first Landau leveEa. 

Turning back to the Gaussian approximation, we note that it has only been 
performed in the bosonized action of the free fermionic term. The interaction part 
is exactly bosonized by just replacing p by B. Considering a two-body scattering 
interaction jpVp, the full bosonized action reads 

(109) 

Note that we were able to transform a fermionic interacting action in an action 
that describes dynamical electric fields and magnetic fluxes. In order to calculate 
current correlation functions we must fix the gauge. In the Coulomb gauge, the 
above action reads 
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47 \^ /Coo ) 4y \k^kJlC^j - k^Ti-IC J 

(110) 

It is now a simple task to evaluate, for example, the density-density correlation 
fmiction: 



k)p{-u;, -k))s = {B{lo, k)B{-u, ~k))s,^^ = 

1 + JCmV{k) 

It is interesting to note, that, in this case, the Gaussian approximation coincides 
with the random-phase approximation (RPA). Of course, equation ( |105| ) is a natural 
framework to systematically improve this approximation. 
The current -current corelation function reads 



{Muj,k)Jk{-uj,-k)) 



S 



eaekm{E^{LJ,k)Ek{-uj,~k))s^,^ = {E^)S^k ~ {E^Ek) 



(P6,, - hk,) {AoAo) + ^^{BB) 
\ / k^ k^ 

{k^5ik - hkk^ 



kiknJCin — fc Tr/C uj kikk JC, 



00 



k^ F fc2 i + /CooT/(fc) 

Finally, we can evaluate the density-current correlation function, obtaining 



{p{iu,k)Jk{-oj,-k))s = €k-^{B{uj,k)Ek{-Lo,-k))s^^^^'^{BB) 

k^ 



(112) 



Lokk /Co 



fc2 l+lCoQV(k) 



(113) 



For the sake of consistency let us check whether current conservation kiJi+ujp = 
holds in the above current correlation functions. Multiplying eq. (|ll2| ) by ki and 



using eq. (113), we obtain 



^2 ^ 

ki{Ji{u!, k)Jk{-uj, -k))s = -r-kk = uj{p{uj, k)Jk{~Lu, -fc))s, (114) 

k^ l + ICoaV{k) 

as it must be. Let us remark that the current conservation in the fermionic version 
is guaranteed thanks to the global U{1) invariance, while in the bosonic counterpart 
it is guaranteed by the topological character of the current. 



5. Conclusions 
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In this paper wc have addressed the question of the universahty of current 
bosonization rules in (2 + l)d non-relativistic fermions with linear and non-linear 
dispersion relations. We have shown that in the linear dispersion relation case, the 
fermionic action bosonizes to a gauge theory. The bosonization process does not 
mix the kinetic and interacting parts of the action, making it possible to analyze 
them individually. We have found that the kinetic part of the fermionic action 
maps into an effective action /rB[j4] that is given by the logarithm of the trans- 
verse Fom-icr transform of the fermionic; determinant. Conversely, the interacting 
part of the action bosonizes by simply replacing the fermionic current by the 
topological bosonic one = e^iupduAp. These rules are universal, in the sense that 
they do not depend on the particular model. Moreover, as the bosonized current 
is topological, it is metric independent and therefore this rule does not depend on 
any local property of the model. 

In the case of a non-linear dispersion relation, we have found that the universal 
rules can be applied if we only consider density-density interactions. The current- 
current interaction destroys universality in this case. The reason for this behavior 
resides in the fact that for non-linear dispersion relations, the fermionic currents 
(for instance ji = not gauge invariant. For this reason, we believe that 

it could be possible to find universal bosonization rules in the general case, provided 
we are able to reformulate the bosonization technique in terms of gauge-invariant 
objects. 

We have applied the functional bosonization formalism to a non-relativistic 
Thirring-like model and have evaluated the spectrum of collective excitations in 
various limits. In the large-mass limit, we were able to exactly calculate this spec- 
trum for arbitrary density-density and current-current interactions. We have found 
collective excitations whose dispersion relation has a gap depending on the partic- 
ular relation between the potentials. An interesting observation is that the parity- 
conserving part of the vacuum polarization tensor (and of course, the expression 
for the collective excitations) is similar to the one- dimensional Tomonaga-Luttinger 
model, up to a coupling constant renormalization. The physical interpretation is 
that the kinetic energy of fermions in two spatial dimensional with a large gap 
(mass m) is frozen. Consequently, the only effect of kinematics is to renormalize 
the coupling constants. Of course, in (2 -|- 1) dimensions there is, in addition to 
the parity-conserving term, a parity-breaking term C2, which has no equivalent in 
(1-1-1) dimension and is the responsible for opening the gap in the excitation spec- 
trum. It is necessary to stress here that these are exact results without any kind 
of approximation other than the consideration of a large gap responsible for the 
insulating face. 

We have also analyzed the massless case in a Gaussian approximation in order 
to evaluate the bosonized action and we have shown that in this approximation 
there are no collective excitations for any density-density potential. However, the 
presence of a current interaction may induce a gapless collective mode with linear 
dispersion relation. 
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Another important application of bosonization is to improve the usual approx- 
imation schemes in many body physics. Since we are bosonizing the interactions 
exactly, all the approximations are made on the bosonized kinetic fermionic action. 
In order to improve the Gaussian approximation, we need a better understanding 
of the structure of the transverse Fourier transform of the fermionic determinant. A 
first step towards this aim was made in refs. 113 and cfl where a topological Chern- 
Simons structure was identified. In this context, the problem of going beyond the 
quadratic approximation in Kb is that of evaluating the Fourier transform of an 
interacting gauge theory. An interesting method to decouple this type interaction 
was recently proposed in ref. El. 

Summarizing, we believe that the existence of universal rules for bosonization 
in two spatial dimensions could help us improve our understanding of strongly cor- 
related systems. In particular, the topological character of the current and the 
topological content of the bosonized action could play a similar role to the one 
that chiral symmetry plays in one dimension, possibly allowing, for instance, the 
calculation of universal properties of two dimensional systems. Thia kind of ap- 
plication is now under development and will be reported elsewhere E3. Of course, 
a definite improvement on the comprehension of the structure of Kb would be to 
discover a fermion-boson mapping between fields (and not between currents) in two 
dimensions. This is one of the most important open questions within this subject. 
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